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The quantum Boltzmann equation (QBE) is the starting point for the calculation of 
the relic density of stable, weakly interacting massive particles (WIMPs) that were once in 
thermal equilibrium. This describes a number of dark matter candidates (such as the lightest 
super symmetric particle, or LSP). A linearized form of the QBE is also used to calculate the 
power spectrum of fluctuations in the cosmic microwave background. We rely on the results 
of these calculations, and expect them to be accurate to a high degree of precision. 

It is therefore somewhat disconcerting to realize that the theoretical derivation of the 
QBE is not at all straightforward, and that there can be important corrections. (For a 
good general discussion of the key issues in the simpler context of nonrelativistic quantum 
mechanics, see [f|j].) Furthermore, I am unaware of any experimental verification of the 
QBE in which the relevant cross sections have been calculated from first principles, with no 
phenomenological input. Thus astrophysics leads us to a new frontier for the QBE. 

In a series of papers, Matsumoto and Yoshimura (MY) have carefully reconsidered the 
derivation of the QBE for WIMP relic densities. They found a surprising result: a previously 
unsuspected correction term that becomes dominant at low temperatures, and changes the 
standard results. Furthermore, this correction can be computed and understood using the 
methods of thermal field theory; the subtleties in the derivation of the QBE do not enter 
into it. 

Simply stated, the result of MY is that the number density n of a gas of WIMPS, each 
with mass M, in thermal equilibrium at a temperature T <C M with a gas of massless 
particles, is given by 

n = (2tt)- 3 / 2 (MT) 3 /V m / t + cA 2 T 6 /M 3 + . . . , (1) 

where A is a coupling constant and c is a numerical factor. The first term represents the 
usual result for noninteracting particles, and carries a Boltzmann suppression factor e~ M ' T \ 
the second term is a loop correction, and it does not have this suppression. Therefore, no 
matter how small A is, at sufficiently low temperatures the second term will dominate. 

Is this result correct? We must first ask for the definition of n. MY use a model of 
a heavy spin-zero boson (represented by a real scalar field (p) interacting with a massless 
spin-zero boson (represented by a real scalar field x)- The scalar potential is 

V{<p, X) = \M\ 2 + ±\^ 4 + ^\X A + Wx 2 - (2) 

We assume that A^ A 2 , A A x , and A x < 1. This hierarchy among the couplings allows 
the light x particles to function as an efficient heat bath for the heavy <p particles. Then, 
for T <C M, MY define the number density n of heavy particles via 



TrH e~ H l T 

T Z- H/T -«Wo> 



(3) 



M x *" M 

where H is the total hamiltonian, and 

H v = \p> 2 + \{Vp) 2 + §MV + counterterms (4) 

is the free-field part of the ip hamiltonian, plus counterterms (some of which involve the x 
field) that are necessary to remove infinities in this composite operator. 
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Eq. ([3D is a highly plausible definition of the heavy particle number density (at low 
temperatures). However, it does not correspond in any obvious way to how this number 
density would be determined experimentally Standard methods all involve a search for 
individual, on-shell ip particles. Real-world examples of this include present-day dark matter 
searches, and measurements of the cosmic microwave background radiation. Therefore, it is 
possible to suspect that the definition (|[) is not appropriate, and this is the reason for the 
surprising lack of Boltzmann suppression in the loop corrections to n. 

In ||, this issue was investigated using a simple Caldeira-Leggett model f|]. In this 
exactly solvable model, a "system" (represented by a harmonic oscillator) is coupled to an 
"environment" (represented by more oscillators) via an "interaction". We found that, at 
low temperatures, the "interaction" energy was always comparable to the "system" energy, 
making the identification of the "system" problematic. 

Consider now a slight variation of the MY model, in which the heavy field is complex, 
and carries a conserved U(l) charge, while the light field remains real and neutral ||. In 
this model, we can study the charge fluctuations in a given volume V. For noninteracting 
particles, Poisson statistics for the number of positively and negatively charged particles in 
V results in 

<Q 2 > = nV. (5) 

We see that the charge fluctuations give us a measurement of the total number of heavy 
particles in a given volume. If Q represents electric charge, we can in principle measure (Q 2 ) 
without tracking individual heavy particles. Furthermore, it seems highly unlikely that weak 
interactions could significantly modify eq. (f|). If (Q 2 ) is either much larger or much smaller 
than nV, then the movements of positive and negative particles into and out of V would have 
to be highly correlated (in order to suppress or enhance the charge fluctuations in V). This 
is inconsistent with the usual notion of a gas of particles that move freely and independently 
between occasional scatterings, and would appear to require strong interactions. 

We therefore define the number density of heavy particles, for T <C M, via eq. (^j). 
We can now compute the 0(X 2 ) loop corrections to (Q 2 ), and see whether or not they are 
Boltzmann suppressed. The answer is: they are. Details of the loop graph evaluation are 
given in ||, but there is a simple argument for the Boltzmann suppression of (Q 2 ) to all 
orders in perturbation theory. Consider the partition function 

Z = Tre^ H ~^ T (6) 

where \i is a chemical potential for the charge Q. We then have 

We will now show that all /i-dependent terms in Z are Boltzmann suppressed. We write 

Z = £ e -(£ Q -MQ.)/T (g) 
a 

where the sum is over a basis of energy and charge eigenstates. States that yield a fi- 
dependent contribution to Z must have Q ^ 0. However, in a weakly coupled theory, the 



(7) 
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lowest-energy state with Q ^ is a state consisting of a single zero-momentum ip parti- 
cle. This is an exact energy eigenstate with energy M; its contribution to Z is obviously 
Boltzmann suppressed. Since all other states with Q ^ have higher energy, their contribu- 
tions to Z are Boltzmann suppressed as well. Eq. (]?]) then implies that (Q 2 ) is Boltzmann 
suppressed. 

In the original model of MY, there is no conserved charge whose fluctuations can be 
measured to determine the number density n of heavy particles. However, is seems likely 
that any definition of n that corresponds to an experimental measurement will yield a 
result that is Boltzmann suppressed. Assuming this, a simple mathematical definition of 
n is n = (l/M)p BS , where p BS represents those terms in the total energy density that have 
a Boltzmann-suppression factor. Any unsuppressed terms in p are to be thought of as 
corrections to the energy density of the gas of light particles. This interpretation is in fact 
the usual one in quantum electrodynamics [|J, where unsuppressed loop corrections to the 
total energy density p (at temperatures well below the electron mass) are attributed to 
self-interactions among photons that arise after integrating the electron-positron field out 
of the functional integral. Recently, this point of view has been investigated in detail by 
Braaten and Jia |7j for the original MY model. They show that, after integrating out the 
heavy (p field to get an effective theory of the light x field alone, the T 6 /M 2 term in p v is 
obtained as the thermal expectation value of a particular higher-dimensional operator in the 
hamiltonian of the effective theory. However, the coefficient of this operator can be changed 
to any value (including zero) by a field redefinition. This shows that no physical significance 
can be attached to the T 6 /M 2 term in p v . 

We therefore conclude that the QBE for WIMP relic densities should not be corrected 
in the manner proposed by MY. However, we should keep in mind that the theoretical 
situation is not as clear as we might like it to be, and so be on the lookout for possible 
further surprises. 
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